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Abstract 

We present ab initio calculations on the effect of in-plane equi-biaxial strain on the structural 
and electronic properties of hypothetical graphene-like GaN monolayer (ML-GaN). It was found 
that ML-GaN got buckled for compressive strain in excess of 7.281 %; buckling parameter increased 
quadratically with compressive strain. The 2D bulk modulus of ML-GaN was found to be smaller 
than that of graphene and graphene-like ML-BN, which reflects weaker bond in ML-GaN. More 
importantly, the band gap and effective masses of charge carriers in ML-GaN were found to be 
tunable by application of in-plane equi-biaxial strain. In particular, when compressive biaxial 
strain of about 3% was reached, a transition from indirect to direct band gap-phase occurred with 
significant change in the value and nature of effective masses of charge carriers; buckling and 
tensile strain reduced the band gap - the band gap reduced to 50 % of its unstrained value at 6.36 
% tensile strain and to 0 eV at an extrapolated tensile strain of 12.72 %, which is well within its 
predicted ultimate tensile strain limit of 16 %. These predictions of strain-engineered electronic 
properties of highly strain sensitive ML-GaN may be exploited in future for potential applica¬ 
tions in strain sensors and other nano-devices such as the nano-electromechanical systems (NEMS). 
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1 Introduction 

Bulk gallium nitride (GaN), which usually 
adopts a wurtzite (WZ) structure, is a wide di¬ 
rect band gap (Eg(T = 0) = 3.51 eV) semicon¬ 
ductor, which makes it suitable for blue lasers 
and light emitting diodes (LEDs) [BIS]. Its al¬ 
loy with indium nitride (InN, Eg(r = 0) = 0.69 
eV), i.e., InjjGai-jjN, has continuously tunable 
band gap that spans the solar energy spectrum 
by varying x mm- Due to its wide and tunable 
band gap, high thermal and mechanical stability, 
GaN is one of the most promising materials for 
applications in photovoltaics, optoelectronic de¬ 
vices, high-temperature microelectronic devices 
and solar cells. Gurrently, band gap engineer¬ 
ing in bulk GaN is mainly done via doping or 
alloying, for instance with In as in In 3 ;Gai_ 3 ;N. 
This method suffers from two fundamental prob¬ 
lems, viz., (I) when In fraction x reaches 0.30, 
the quantum efficiency of LEDs drops signifi¬ 
cantly to undesired levels; (2) when x reaches 
0.32, a phase separation occurs because of the 
large atomic size mismatch of Ga and In [3l[3]. 
To overcome these problems, alternative meth¬ 
ods of band gap tuning in GaN through equib- 
iaxial in-plain strains have been proposed re¬ 
cently [3j as a method of ’’strain engineering”. 
In the context of GaN based nano-materials, 
nanowires, nanotubes, and nanospirals of GaN 
have been synthesized which showed great po¬ 
tential for fabricating wide-spectrum LEDs and 
other nano-scale devices m\- In particular, a 
density functional theory (DFT) study in 2006 
[8] predicted that when the layer number of 
(OOOI)-oriented WZ materials (e.g., AIN, BeO, 
GaN, SiC, ZnO, and ZnS) is small, the wurtzite 
structures transform into a new form of stable 
hexagonal graphite-like or hexagonal boron ni¬ 
tride (BN)-hke structure. This prediction has 


been confirmed by experiments in respect of ZnO 
in 2007 [9]. Recently the thickness range of 
graphitic WZ films (e.g., AIN, GaN, InN, SiG, 
BeO, ZnO) are theoretically shown to depend 
on strain and can be extended to much thicker 
films by epitaxial strain [TO]. In particular for 
GaN film, it has been shown [TO] that (i) a 10- 
atomic-layer (000I)/(000I) GaN film transforms 
into a more stable 5 monolayer (ML) graphitic 
film having no destabilizing dipole moments in 
its surface, (ii) For GaN films beyond lO-layers, 
the WZ structure tends to be stable due to the 
charge transfer from the (0001) bilayer to the 
(0001) bilayer, counteracting the intrinsic dipole 
perpendicular to the film, (hi) with 10% biaxial 
tensile strain, the graphitic phase of GaN was 
found to be stable up to 24 layers and (iv) the 
band gap of the stable strained graphitic films 
are also shown to be tunable over a wide range 
either above or below that of the bulk WZ phase. 

The study of 2D crystals is an emerging field 
of research inspired by the recent phenomenal 
growth in the research on graphene (a one atom- 
thick nanocrystal of carbon atoms tightly-bound 
by the sp^-hybridized bonds in a 2D hexago¬ 
nal lattice) which promises many novel applica¬ 
tions mHni. Representative samples of some 
2D/quasi-2D nanocrystals that have been syn¬ 
thesized recently include ZnO [9], BN [TTtili] . 

M 0 S 2 [IIHIll, MoSe 2 [nHHI, BDTes [IIHIl], 

Si |I5j . With rapid advancement of synthe¬ 
sis techniques, the synthesis of GaN monolayer 
(ML-GaN) is expected in near future. There¬ 
fore, understanding and controlling the elec¬ 
tronic properties of ML-GaN is of much signif¬ 
icance both from fundamental and application 
point of view. Although several theoretical stud¬ 
ies have been reported |I6It2Ij to this end, the ef¬ 
fect of strain on the structural and the electronic 
properties of ML-GaN remains unclear. Since 
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Figure 1: Buckling probe of planar (PL) ML- 
GaN based on the energy minimization proce¬ 
dure with two atoms per unit supercell. Insets 
show the ball-stick model for (a) the side view 
of buckled (BL) ML-GaN, (b) the side view of 
planar ML-GaN and (c) the top-down view of 
both PL and BL ML-GaN. In BL ML-GaN, Ga 
atoms (large balls, red) and N atoms (small balls, 
green) are in two different parallel planes; buck¬ 
ling parameter is the perpendicular distance be¬ 
tween those parallel planes and A = OA for PL 
ML-GaN. 


strain-engineering is an important method to tai¬ 
lor the structural and electronic properties of 
nanomaterials |22H24j , using DFT we have inves¬ 
tigated the effect of equi-biaxial inplane strain on 
ML-GaN. Here we report the results our study 
which simulates an ideal experimental situation 
in which ML-GaN is supported on a flat stretch¬ 
able substrate. 


The calculations have been performed by us¬ 
ing the density functional theory (DFT) based 
full-potential (linearized) augmented plane wave 
plus local orbital (FP-(L)APW-l-lo) method [25] . 
which is a descendant of FP-LAPW method [26]. 
We use the elk-code m and the Perdew-Zunger 
variant of local density approximation (LDA) 
[28] . the accuracy of which have been success¬ 
fully validated in our previous studies [291 - I36] 
of graphene and some graphene-like 2D nano¬ 
crystals. For plane-wave expansion in the in¬ 
terstitial region, we have chosen |G -|- k|niax x 
Rmt = 8.5 (Rmt is the smallest mufRn-tin radius 
in the unit cell) for deciding the plane-wave cut¬ 
off. The Monkhorst-Pack I37]k -point grid size of 
30x30x 1 was chosen for all calculations. The to¬ 
tal energy was converged within 2//eV/atom. We 
simulate the 2D-hexagonal structure of ML-GaN 
as a 3D-hexagonal supercell with a large value of 
c-parameter (= |c| = 40 a.u.) and two atoms per 
unit supercell. The application of in-plane equi- 
biaxial strains about ±10% was simulated by 
varying the in-plane lattice parameter a(= |a| = 
|6|); e = [(a—ao)/ao] x 100, where oq is the ground 
state in-plane lattice constant of ML-GaN. Our 
consideration of such large values (±10%) of 
strain is inspired by (i) the first principles study 
of mechanical properties of ML-GaN [2T], where 
ML-GaN is shown to have an ultimate strain 
value of 0.16 under biaxial deformation, (ii) the 
recent theoretical studies on strained GaN films, 
where strains in the range of (—5% to±10%) [lOj 
and (—20% to ± 20%) |1| have been considered 
and (iii) the fact that graphene (whose honey¬ 
comb structure the ML-GaN adopts) is experi¬ 
mentally demonstrated |3S| to sustain in-plane 
tensile elastic strain in excess of 20%, although 
we do not yet clearly comprehend the maximum 
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compressive strain graphene can sustain while 
preserving its planar structure. However, for the 
sake of symmetry in the compression variation, 
and as a matter of curiosity, we have considered 
compressive strains up to 10%. For the struc¬ 
tural stability of ML-GaN, we have considered 
two possible structures of ML-GaN (with two 
atoms per unit cell) as presented in Figure 1. 

3 Results and Discussions 

3.1 Structural and Mechanical Prop¬ 
erties 

For an assumed flat ML-GaN, our calculated 
LDA value of oq = 3.1561 A corresponds to 
the Ga-N bond length = ao/\/3 = 1.822 

A, which is close to the reported LDA value of 
^Ga-N ~ ^ §ahin et al. [16], GGA val¬ 

ues of 1.87 A by Chen et al. [T7], 1.83 A by 
Xiao et al. [18], 1.852 A by Peng et al. [21]. It 
is well known that usually LDA underestimates 
and GGA overestimates the lattice constant. In 
our previous studies of graphene and silicene [29] 
and germanene [3D], we have demonstrated that 
the value of c—parameter chosen in the construc¬ 
tion of supercell for simulation of 2D hexagonal 
structures also affects the value of in-plane lat¬ 
tice constant oq: larger value of c—parameter 
yields a smaller value of oq. Since we use a differ¬ 
ent method and a different value of c—parameter, 
the slight disagreement of our result on oq with 
other theoretical results is well understood and 
acceptable. Our assumption on the flat ML-GaN 
structure was tested correct by the calculated re¬ 
sults depicted in Figure 1, based on the princi¬ 
ple of minimum energy for the stable structure. 
This corroborates the results of previous theo¬ 
retical studies [HHinittSHiaiii] based on differ¬ 
ent methods and the stability of flat ML-GaN is 


attributed to the strong in-plane sp^ hybridized 
bonds between the Ga and N atoms. 

The calculated 2D bulk modulus (also called 
in-plane Young’s modulus) of ML-GaN, B 2 Y) = 
A{d'^E/dA^)\A^i„ = 109.45 N/m, (where A is 
the area of the periodic cell of the 2D lat¬ 
tice and Amin is the area with minimum en¬ 
ergy) is in good agreement with reported results: 
110 N/m [16], 109.4 N/m [2T]. This value is 
smaller than the corresponding calculated values 
of graphene (214.41 N/m) and ML-BN (181.91 
N/m) [39], reflecting the weaker bonds in ML- 
GaN than those in graphene and ML-BN. There¬ 
fore, ML-GaN is expected to be useful in ap¬ 
plications requiring materials with less tougher 
behavior than graphene and ML-BN. Further, 
since i? 2 p) (ML-GaN) « (Graphene) as 

seen above, we expect ML-GaN to sustain elastic 
tensile strains ~ 10% considering the fact that 
graphene sustains elastic tensile strains in excess 
of 20% |38j . However, to know the compressive 
strain value that ML-GaN can sustain without 
change in its planar structure, we have made 
some study by probing the buckling in ML-GaN 
at different values of compressive strains. The 
results of this probe for four values of compres¬ 
sive strains, viz., e = —7%, —8%, —9% and—10%, 
are shown in Figure 2, in which it is seen that 
ML-GaN sustains compressive strain of 7% with¬ 
out buckling and with higher compressive strain 
values of 8%, 9% and 10%, buckling parameter 
assumes nonzero and higher values. From these 
data we inferred that there is a critical compres¬ 
sive strain value beyond which buckling starts. 
To determine this critical point, we have plotted 
the buckling parameters calculated at six com¬ 
pressive strains of 7.5%, 8.0%, 8.5%, 9.0%, 9.5% 
and 10% in Figure 3, where our data fit well with 
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Figure 3: Variation of buckling parameter of 
ML-GaN with absolute value of compressive 
strain. 
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Figure 2: Buckling probe of ML-GaN at four 
different values of compressive strain of (a) 7%, 
(b) 8%, (c) 9% and (d) 10%. 


the following second order polynomial, viz., 

A(|e|) = -3.42492-b0.6961|e| -0.03104|ep (1) 

where A is in A and |e| represents the absolute 
value of compressive strain in %. One of the solu¬ 
tions of Eq.(l) for A(|e|) = 0 A is |e| = 7.281% 
which is the critical compressive strain beyond 
which ML-GaN becomes buckled. This is in 
agreement with our calculated result of Figure 
2(a), where we found no buckling at compressive 
strain of 7%. Eq.(l), which is valid for compres¬ 
sive strain value |e| < 10%, may be useful for 
calculating the buckling parameter correspond¬ 
ing to a particular value of compressive strain 
|e| < 10%, which in turn may be used to calcu¬ 
late the electronic properties of ML-GaN under 
compressive strains 10% < |e| < 7.281%. So 
we have used this approach to calculate the elec¬ 
tronic properties of buckled ML-GaN under high 
compressive strain values |e| < 10%, which we 
present in the next section. 
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Figure 4; Bands and total DOS of unstrained 
ML-GaN within LDA. is the Fermi energy. 


3.2 Electronic Properties 

The electronic band structures and total den¬ 
sity of states (TDOS) plots of unstrained ML- 
GaN are depicted in Figure 4. As seen in Fig¬ 
ure 4, unstrained ML-GaN is an indirect band 
gap (Eg = 2.462 eV {K —>• F), LDA value) 
semiconductor with conduction band minimum 
(GBM) located at the F point and valence band 
maximum (VBM) located at the K point of the 
hexagonal Brillouin Zone (BZ). However, the ac¬ 
tual band gap is surely larger than this since 
LDA is well known to underestimate the gap. 
Our calculated LDA band gap of 2.462 eV may 
be compared with reported LDA value of 2.27 
eV by §ahin et al. [16], GGA values of 1.95 A 
by Ghen et al. 2.46 eV by Xiao et al. m- 
Using the GW approximations, the GW band 
gap of ML-GaN have been estimated at 5.00 eV 
(LDA-l-GWo) by §ahin et al. [TB] and 4.14 eV 
(GGA-I-GW) by Ghen et al. |T7|, which showed 
that ML-GaN is a wide indirect band gap semi¬ 
conductor. Although both LDA and GGA do 
not yield the band gap correctly, these are power¬ 
ful enough to predict a correct trend in variation 
of the gap dlttZlllolllI]. Since we focus on the 
trend as well as the relative variation of Eg rather 
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Eigure 5: Relative variation in band gap 
6Eg (with respect to the unstrained band gap 
Eg(0) = 2.462 eV) with in-plane homogeneous 
biaxial strain e. 


than its absolute value, we employed the compu¬ 
tationally simpler and less time-consuming LDA 
for this study. 

We simulated the application of in-plane equi- 
biaxial strain (i.e., hydrostatic strain) to ML- 
GaN by varying the in-plane lattice parameter a 
from 2.84 A to 3.46 A with small increments of 
0.02 A. With our calculated value of oq = 3.1561 
A for unstrained ML-GaN, we thus simulated 
the strains in the range (—10% to 4-9.63%). Our 
calculated results on the in-plane equi-biaxial 
strain-induced modifications of the band gap of 
ML-GaN are depicted in Eigure 5, which shows 
both the trend as well as the relative change in 
Eg defined by 

^Eg = [(Eg(e) - Eg(0))/Eg(0)] X 100 (%) (2) 

with respect to the unstrained band gap Eg(0) = 
2.462 eV (the relative change is expected to be 
more meaningful and free from computational 
method dependent errors). For 3.08 A< a < 
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3.46 A (-3.045% < e < 9.629%), Eg remains 
indirect {K —)• E) and its variation with strain 
is linear; for 2.84 A< a < 3.06 A (—10% < 
e < 3.045%), Eg remains direct (E —)• E) and its 
variation is non-linear with strain in the planar 
phase and almost linearly decreasing in the buck¬ 
led phase. The band gap variation is asymmetric 
with respect to e not only in the trend but also in 
the character (transition from indirect to direct 
band gap). For tensile strains the relative change 
in Eg is very much pronounced; Eg dramatically 
drops to 50% ((5Eg = —50%) of its unstrained 
value at e = 6.36% (obtained by linear interpola¬ 
tion) and to 0 eV (5Eg = —100%) at e = 12.72% 
(obtained by linear extrapolation). In contrast, 
for high compressive strain, the reduction in Eg 
is about 44.78% at e = —10%. More interest¬ 
ingly, a transition from indirect to direct band 
gap occurs at about 3% compressive strain as 
noted in Figure 5 for Eg(a = 3.08A) = 2.863 
eV {K —>• E) and Eg(a = 3.06A) = 2.943 eV 
(E —)• E). Further, the direct band gap (which 
is larger than pristine indirect band gap) non- 
linearly saturates ((5Eg = -|-22%) with increasing 
compressive strain until ML-GaN assumes buck¬ 
ling, which in turn induces a reduction in Eg 
as seen in Figure 5. This property of ML-GaN 
may be exploited for non-linear high-frequency 
optical applications in future. It is to be noted 
that, except for the occurrence of the direct to 
indirect gap-phase transition at about 3% com¬ 
pressive strain and the positions of GBM and 
VBM, the present trend in the variation of Eg 
with e is similar to those recently predicted for 
bi-axially strained bilayer GaN (2-ML-GaN) [lOj 
and ultrathin GaN films in graphitic phase; 
the positions of VBM and CBM in 2-ML-GaN 
are not reported in ref. m, but in graphitic GaN 
films, the VBM is at H point and CBM is posi¬ 
tioned at E point of the 3D hexagonal BZ [4]. 


Thus the electronic properties of ML-GaN, are 
different from those of few layer graphitic GaN 
films. This is not surprising, as we know that 
the electronic properties of monolayer, bilayer 
and few-layer graphene are all different. A re¬ 
lated theoretical study [32] of bi-axially strained 
direct band gap ML-BN (corroborated by our 
recent study m). showed strain induced gap- 
phase transition from direct to indirect band gap 
at e = —1.53%, and Eg follows linear variation 
for strains below and above e = —1.53%, which 
substantiate our confidence on the present re¬ 
sults. Since the band gap of ML-GaN is very 
sensitive to strain, it is a potential candidate for 
new-generation strain sensor. 

The nature of and trend in the strain-induced 
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direct gap-phase, calculated at the band edges 
at the special points K and F as appropriate 
for the two gap-phases, are shown in Figure 6. 
For unstrained ML-GaN, we calculated me.(T) = 
0.266mo and mh{K) = 1.157mo (mo is electron’s 
rest mass), which are larger than the band edge 
electron effective mass of 0.20mo for the bulk 
wurtzite GaN mm- As seen in Figure 6, un¬ 
like the case with electrons, the character and 
the value of effective masses of holes abruptly 
change at about e = —3%. The discontinuity 
in the variation of the effective masses of holes 
is due the change of the gap-phase which makes 
the curvature of the VBM at K point in the indi¬ 
rect gap-phase much smaller than the curvature 
of VBM at the F point in the direct gap-phase. 
Since the mobility of electrons and holes depend 
on their effective masses, the carrier mobility in 
ML-GaN can also be controlled by strain engi¬ 
neering. 

4 Conclusion 

In summary, we have theoretically explored 
the effects of in-plane equi-biaxial strains on 
the structural and electronic properties of 
ML-GaN in graphene-like planar honeycomb 
structure employing DFT based all electron 
full-potential calculations. In this study we 
found that ML-GaN gets buckled for compres¬ 
sive strain in excess of 7.281 %. Within our 
considered compressive strain limit of 10%, 
buckling parameter varies quadratically with 
compressive strain with an increasing trend 
with increasing compression. The calculated 
value of 2D bulk modulus of ML-GaN reflects 
weaker bond in ML-GaN in comparison to that 
in graphene and ML-BN. More importantly, 
our calculations show that the band gap and 


effective masses of charge carriers in ML-GaN 
can be significantly tuned and manipulated by 
application of in-plane equi-biaxial strain. In 
particular, when compressive biaxial strain of 
about 3% is reached, a transition from indirect 
to direct band gap-phase is predicted to occur 
with significant change in the value and nature 
of effective masses of charge carriers; the band 
gap reduces to 50% of its unstrained value 
at 6.36% tensile strain and to 0 eV at an 
extrapolated tensile strain of 12.72%. These 
interesting predictions of strain-engineered 
electronic properties of highly strain sensitive 
ML-GaN may be exploited in future for novel 
device applications such as nano-mechatronics, 
strain sensor, nano-electromechanical systems 
(NEMS) and nano-opto-mechanical systems 
(NOMS). 
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